Path Integral Approach to Strongly Nonlinear Composites 
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We study strongly nonlinear disordered media using a functional method. We solve exactly the 
problem of a nonlinear impurity in a linear host and we obtain a Bruggeman-like formula for the 
effective nonlinear susceptibility. This formula reduces to the usual Bruggeman effective medium 
approximation in the linear case and has the following features: (i) It reproduces the weak contrast 
expansion to the second order and (ii) the effective medium exponent near the percolation threshold 
are s = 1, t = 1 + n, where k is the nonlinearity exponent. Finally, we give analytical expressions 
for previously numerically calculated quantities. 



The study of the properties of linear heterogeneous me- 
dia (composites, suspensions) has been the subject of an 
intense activity for already fifty years (see the reviews 
and H ) . More recently there has been a great interest in 
non-linear media The nonlinear composites are im- 
portant for technology and also from a fundamental point 
of view, for which it is important and challenging to un- 
derstand the interplay between nonlinearity and disorder. 
There are essentially two types of nonlinear media, (a) 
Weak nonlinearity: the nonlinearity is small compared to 
the linear term. This case was studied by many authors 
and is now relatively well understood j^-|ll]].(b) Strong 
nonlinearity: the nonlinearity is here the dominant term 
and this can happen in essentially two situations. First, 
there can be a sharp threshold between two different be- 
haviors and this case can model the phenomena of frac- 
ture or dielectric breakdown |l2|]lg| l. Second, the consti- 
tutive law can be a pure power law of the form 



3 = X\E\ K E 



(1) 



where j is the current and E the electric field. This be- 
havior can be observed in a dielectric illuminated by a 
Laser [Q] when the multiphoton processes dominates and 
the usual linear approximation completely breaks down. 
Certain cermets resistors, ZnO based varistors 14,191 or 
disordered alloys |]l6| can also display this behavior (|1|). 

In a disordered medium, the nonlinear susceptibilities 
X can fluctuate from point to point and one is interested 
in the macroscopic effective behavior of such a medium. 
If the nonlinearity exponent k is the same for all phases 
of the medium, then the effective nonlinear susceptibil- 
ity is well-defined and is given by jo = Xe\Eo\ K Eo where 
jo and Eo are the macroscopic current and electric field 
respectively. It is difficult to evaluate Xe and devising a 
reliable method to compute Xe would allow one to study 
a variety of other problems such as fracture or dielectric 
breakdown. 

In this strongly nonlinear case, Blumenfeld and 
Bergman obtained the weak contrast expansion to sec- 
ond order ]l7| . This expansion was recovered by means 
of a path integral method 18 1. The dilute limit was stud- 
ied in |19J20tl. Problems arise when one tries to find an 



effective medium approximation (EMA) for this type of 
media. A good EMA should satisfy the two following 
criteria, (i) It should reproduces the weak-contrast ex- 
pansion (at least up to the second order) and the dilute 
limit (although this last condition is probably very dif- 
ficult to fullfill for strongly nonlinear media), (ii) Close 
to the percolation threshold p Cl the effective nonlinear 
susceptibility Xe is described by two exponents: For a 
metal/insulator mixture, one has pl|,p2[ 



Xe ~ (P - Pc) 



t(«) 



(2) 



where p is the proportion of the conducting component. 
For a superconductor/metal mixture, one expects 



Xe ~ (Pc - P) 



(3) 



where p is here the proportion of the superconduct- 
ing component. Diffcrents values for these exponents 
were proposed. In p2fl , the effective medium values are 
t(n) = 1 + k and s(k) = 1, and in p3| , p4fl , the exponents 
are t(n) = s(n) = 1 + k/2. In both cases, the crossover 
exponent <f> — s + 1 is equal to 2 + k. These two sets 
of exponents satisfy the duality relation for d — 2 [5l| : 
t(j) = 75(1/7) where 7 = 1+k. So far, numerical results 
p5| and series analysis [§6| suggests that for d = 2 the 
exponents s and t are different, ruling out s — t = 1 + k/2 
although further numerical studies are necessary to make 
a definitive statement. An acceptable EMA should pre- 
dict such kind of values. 

We can distinguish two differents classes of approaches 
to this problem. A first approach |t|,|24],^7] |u| consists in 
expressing the effective nonlinear susceptibility in terms 
of the averaged electric field in each component. In 
p4] , ^7| -|2g[| , a kind of a "decoupling approximation" is 
proposed for calculating these fields, and one obtains 
a set of coupled equations which is solved numerically. 
Although the agreement with numerical simulations is 
generally fairly good, there are a few drawbacks to this 
method. In particular, this method relies quite heavily 
on numerics and it is difficult to check some analytical 
properties. Moreover, the weak-contrast expansion (con- 
dition (i)) is usually not recovered and the exponents are 



1 



difficult to estimate. In particular, the mean-field theory 
proposed in ]24| does not reproduce the weak contrast ex- 
pansion but instead the lower bound established by Ponte 
Castaheda et al fij] . In this case [M , the values of the ex- 
ponents are s — t — n/2 + 1. In another series of papers 
po| , ^l| , the nonlinear host is linearized up to the second 
order and the local electric fields are computed in a self- 
consistent way. With this method, Ponte Castaneda and 
Kailasam [ [30| proposed an effective medium approxima- 
tion which reproduces the weak contrast expansion, but 
for which the exponents are difficult to estimate. 

The second approach is in the Bruggeman spirit |3^ | 
and consists essentially in considering an impurity in an 
effective host. Bruggeman's theory was reformulated in 
order to apply to this problem p2fl an d it was further in- 
vestigated by different authors ||25| , |33| -|35f . This approach 
predicts the effective medium values s — 1 and t = 1 + k 
and reproduces the weak contrast expansion to second 
order. However, all the studies so far are mostly numer- 
ical and we propose here the analytical solution. The 
obtained formula is relatively simple and satisfies con- 
ditions (i) and (ii). We give the analytical expressions 
for numerically estimated quantities |2qj33| , |35|| . The only 
drawback of our result is that the percolation threshold 
depends on k (in the same way as in |35|]). In the discus- 
sion, we address this point and propose a possible way to 
correct this wrong behavior. 

The constitutive relation is j = x(r) K and the 
local energy density associated to it is 



w(r, E(r)) 



X(r) 
k + 2 1 



E 



K + 2 



(4) 



E is the applied field and for a heterogeneous medium 
the quantitie x( r ) a ^ point r is distributed according to a 
binary law (a generalization of our method to other types 
of disorder should be without problems) 



p(x = x(r)) = pS{x - xi) + q$(x - X2) 



(5) 



The total dissipated energy is given by W* — 
Xe|£'o| K+2 / ( K + 2), and can be expressed as a constrained 
minimum 



W* = 



l E=E ,E= 



-V0 



d d rw{r,E{r))) (6) 



Here, we have assumed that W* is a self-averaging quan- 
tity in the thermodynamic limit, which allows us to com- 
pute the average over the disorder (the brackets (•) de- 
note the average over disorder or equivalently, the spatial 
average). The minimum in Eq. (^) can be written with 
the help of path integrals 111] 



The important quantity to study is thus the 'partition 
function' 



Z = / VEe 



-/3H 



(8) 



In a preceeding paper |18| , we made a perturbation 
expansion up to the second order in disorder and we re- 
covered known results We also showed in another 
paper ]3r| how to recover Bruggeman's approximation in 
the functional framework and we recall briefly the idea. 
We start from the expression (p|) and we add and sub- 
stract a Gaussian ansatz 7Yq = J d d rwo(E(r)) 



■f)]7; e -l3'Ko e -l3('H-'Ha) 



(9) 



We expand the second exponential and resum it keeping 
only the contribution at the same point 



-p(H-H ) _ 



E 

fc=l 



d d r(w(r, E{r)) - w {E(r))) 



/OO 
d d r l w (r, E(r)) - w (E{r))} k 

k=l 



^ re -/9[^(r,E(r))-iuo(J5(r))] 



(10) 



(here and in the following, we omit unimportant volume 
factors and cut-offs). The partition function Z is thus 
given by 



d d r 



-f) Ee -f3Ho e -(3[w(r^E(r))~w Q (E(r))] 



(11) 



The following physical picture can be associated with 
this approximation. The background is described by 
7io and at point r there is an impurity described by 
w — wq. The ideal case would be to take a nonlinear 
background described by an effective nonlinear suscepti- 
bility Ho = / |-E| K+2 and a nonlinear impurity, which 
is so far impossible to compute. We thus have to resort 
to a further approximation. The averaged value of the 
electric field is fixed and given by E — Eq. It is thus 
reasonable to assume that the electric field in the back- 
ground will not fluctuate too much (at least far from the 
impurity) and we can expand the nonlinear background 
around E{r) = E . We write E(r) = E + e{r) and 
expand up to the second order in e 



Wq 



Xe 



E£ +z + E*E ■ e(r) 
EoiEoj . 



i^^ ei (r)(^+«=^) ej (r) 



(12) 



W* = lim 



P- 



- - I VE, ■ 

P 



(7) 



where the "Hamiltonian" is H = J d d rw(r, E(r)) and 
where the measure is VE = V{E, (f>)5(E- E )5(E+V(f>). 



while we keep the exact expression for w — \E\ K+2 . 
The final picture is then the following: we compute ex- 
actly the perturbation induced by a nonlinear impurity 
in the nonlinear linearized effective medium. In addi- 
tion to allow calculations, this scheme ensures that the 



2 



weak-contrast expansion will be recovered. The numeri- 
cal study of this problem can be found in [^5 33 35[|. 

The path integral ( |Tl| ) together with the approxima- 
tion raj) can be computed and after some calculations, 
one is led to 



(13) 



with 



M (u) = <^±-u K+2 



■ XeE^U ■ E 



XeES -(Eo + u).(l- 1 1 + K ^^).(E 0+ u) (14) 



Ex 



where u is a d-dimensional vector. The quantity I is 
given by 



I(d, k) 



Sd~ 



d6sin d - 2 9 



cos 



1 + K COS 2 6 



(15) 



where Sd — 2n d /T(d/2) is the surface of the d- 
dimensional sphere. In the linear case, I(d, k = 0) = 1/d, 
and for d = 2, 3, one has p7|Ji^] 



1(2, k) 
1(3, k) 



1 



1 



(16) 
(17) 



We are interested in the large (3 limit [Eq. ([?])] so we can 
apply the saddle-point method to the integral (|l3|). The 
saddle point u* is parallel to Eq and is given by 



-S (k- i) 



1 



(18) 



where here x = Xe/X- One can note that u* is the elec- 
tric field in the nonlinear impurity embedded in the lin- 
earized effective homogeneous host. The effective energy 
W* = -4 (In Z) is then given by W* = W + AW where 

Wq = -^E£ +2 . The natural self-consistent condition 
AW — can be rewritten for the binary disorder (0) 
under the form [E3.p5f 



Xi X2 







(19) 



where 



1 V 



K + 2 



X K + 2 
1 



— 2 + i + M*) 

2[1 + v*(x)] 2 (1-~ + k) 



(20) 



where = u*(x)/Eq is given by Eq. (|18|). This equa- 

tion ( |l9| ) [together with ( fl8| ) and (pc|)] is our main result, 
and we will now discuss it. 



In the linear case (k = 0), one can easily check that 
( |l9| ) reduces to Bruggeman's Q equation. Moreover, in 
the one-dimensional case, one recovers the exact result 

1/Xe = (1/X^) K+1 . 

As expected, the weak contrast expansion is recovered 
up to the second order, namely 



Xe 



(x>- 



2(X) 



(6 X 2 )I(d, K ) 



(21) 



This fact is not surprising since we used as an ansatz Ho 
the nonlinear effective medium linearized up to second 
order. 

Our approximation will not reproduce the exact dilute 
limit, since the nonlinear host is linearized. Instead, we 
will obtain the following expansion to the first order in 
concentration (exact for k = 0) 



Xi 



Xe ~Xi+<7(K + 2)xi/(— ) 



X2 



(22) 



where q is the fraction of component \2- 

The critical behavior is determined by f(x) for x ~ 
and x —>■ oo J2^|35|]. One obtains from Eqs. (18 2(j) 
f(x ~ 0) ~ /(0) - a(n + l)x~ where 



/(0) 



1 



1 

21 



K 

2 



-( — + «) 
2 V / ^ 



■■ + 2 

■■ + 1 



(23) 



For a; — > oo, one has /(a;) ~ /(oo) + - where 



/(oo) 



6 = 



2^1-1 



1 + 2(-) + k) 
2(1 - } + «) 



(24) 
(25) 



For k = 0, one recovers the known exact expressions 
a = d 2 , b = [d/(d - l)] 2 , /(0) = d and 
/(oo) = — d/[2(d — 1)] (the function / is defined up to 
a constant factor, and there is a global additional factor 
1/(k + 2) in our result). These four different coefficients 
were estimated numerically [^,^,^5) . 
The percolation threshold is given by 



Pc 



/(oo) 



1 



/(oo)-/(0) K + 2(-\ + K) 



(26) 



and the exponents are t(n) = 1 + n and s(k) = 1. In 
Fig. 1, we compare for d — 3 this exact expression for 
p c to numerical results |^o] , ^3|,^5| . It thus seems that the 
variational method used in ]25| , |33| does not lead to the 
correct values of the electric field around the impurity. 
The fact that the percolation threshold depends on k is 
the bad feature of this approximation. However, for k 
not too large, or for a contrast not too high this approx- 
imation works well. When one adds a linear background 
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of effective conductivity a e , an additional dimcnsionlcss 
factor A = a e /xeE$ is introduced in the equations. Close 
to p c , the effective conductivity behaves as a e ~ Ap 1 and 
the effective nonlinear susceptibility as Xe ~ Ap 1+K . The 
factor A is then diverging and the behavior of / is mod- 
ified. One then recovers the Bruggeman value for the 
percolation threshold p c = 1 jd. ft thus seems that in the 
one-impurity scheme a linear background is necessary to 
"regularize" the wrong behavior of p c . It might be a way 
to obtain an EMA satisfying conditions (i) and (ii), and 
which gives a correct value for p c . 
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FIG. 1. p c versus k for d = 3. The line is the analytical 
expression (26), the circles and the diamonds represent the 
numerical results from [35] and [25,33] respectively. 
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